ID spirals: is multi stability essential? 
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Abstract 

The origin of ID spirals or antisymmetric ID pulses is explaind so far on the basis 
of multistability of spatially inhomogeneous and temporally oscillatory phases and 
so called nonvariational effects. Thus, coupled amplitude equations which are valid 
near a co-dimention 2 point and provides with the necessary multistable environment 
are commonly in use in the numerical calculations to generate such structures. In the 
present work we analytically show that a complex Ginzburg-Landau type amplitude 
equation which is valid in the Hopf region of phase space near an instability threshold 
admits solutions like antisymmetric pulses traveling in alternate directions from a 
core. The pulses can have well defined spatial profile like Hermite polynomial of 
order unity on the side of the direction of its motion. Thus, a pulse moving in 
the positive direction of x-axis from near a core at origin will be a peak while that 
moving in the negative direction is a trough and a temporal oscillation of such 
structures would give the impression of alternative peak and troughs are moving 
appart being generated from a common core.. 
PACS number (s): 87.10.+e, 47.70.Fw 

One dimensional spiral pattern, experimentally observed in CIMA reaction is a fasci- 
nating phenomenon of nonlinear waves. As described in the experimental part of pQ that 
at first a transition between a stationary periodic and a traveling wave like Hopf mode 
is observed by lowering the starch concentration in the chemical reactor. Now, keeping 
a low starch concentration and by increasing the malonic acid concentration a similar 
transition takes place which makes an almost ID Turing state loose its stability to a Hopf 
state when malonic acid concentration is doubled. It has been observed in this experiment 
that when the Hopf state takes up from the Turing state very often there remain a few 
spots, reminiscent of the previous Turing state, acting as the source of one dimensional 
anti synchronous wave trains. Bands of maximum intensity spreads alternatively toward 
right and left directions with a time delay. As experiment suggests, this phenomenon is 
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essentially one dimensional and variation of two parameters like concentrations of starch 
and malonic acid indicates proximity to a co-dimension 2 point. On the basis of this 
observation a theory has been developed being based on multisability and non variational 
effects. Such localized structures near a Hopf- Turing instability boundary have been ana- 
lyzed in many subsequent papers [2-9] . A common belief is that in the range of parameters 
far away from depinning transitions a pinned Turing state in a global traveling wave phase 
helps sustain ID spirals by accepting one on one side of it while leaving the other on the 
other side. In such a context it is interesting to see if there are other mechanism that can 
explain the generation of antisymmetric ID waves from a core without the consideration 
of multi stability of states. 

In the present paper we are going to show an alternative scenario where the multi 
stability is not an essential ingredient to make asymmetric outgoing pulses generate from 
a source. Here, such ID asymmetric moving pulses are shown as a steady state solution 
on uniformly moving frame to a complex Ginzburg-Landau type amplitude equation. In 
a purely Hopf environment where such an equation is applicable and there is no coupling 
to typical Turing type amplitude equation, small sources can exist which separates two 
regions where alternative peak and trough of chemical concentrations can move in outward 
directions from the source. In the present work we emphasize on the fact that generation of 
ID spirals can definitely be understood on the basis of a single amplitude equation analysis 
which is generally ignored with a priori consideration of such things as non variational. 
On a moving frame of reference the complex Ginzburg-Landau type amplitude equation, 
applicable to pure Hopf regime, admits steady local structures which have the spatial 
profile like quantum mechanical harmonic oscillator steady states (Hermite polynomials). 
To get such linear solutions the nonlinearity is considered negligible which in tern allows 
only two orders of the polynomial namely zeroth and the first. These two orders are quite 
different in symmetry and helps explain the antisymmetric effect. These local structures 
can oscillate when formed in a Hopf environment. It would be shown that on the moving 
frame the origins of such localized structures are a little shifted on a direction toward 
which it moves where the amount of shift is proportional to the velocity of the frame. 
Thus, by considering two moving frames one toward right and the other toward left of 
some central core one can easily think about a pair of such localized structures moving 
opposite to each other starting from two points a little right and a little left of the actual 
origin on lab frame. Moreover, it would be shown that such structures are well defined 
on the upper (lower) side of their center of a right (left) words moving frame. The inner 
sides of both of them being less well defined we would better observe the positive x part 
of the rightwards one and the negative x part of the left moving one and this can give 
the impression of an antisymmetric wave being generated from a core, thus, when a 
concentration peak is moving toward right a trough moves toward left with the same 
velocity. 

Let us take the amplitude equation in the form 



dH 
~dt 



(fj, r + ijJ,i)H + (D r + iDi) 



d 2 H 



((3 r + tP t )\H\ 2 H 



dx 2 



where H is the slow and large scale dependent amplitude of the Hopf mode. 
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Now, Taking a Hopf phase of the form H = He tuJt where H is a function of space and 
time we get, 

dH d 2 H 

iuH+ — = {fi r + t f i l )H + (Dr + iDi)—^ - (ft + iPi)\H\ 2 H 

After shifting to moving frame with the transformation y = x — vt where v is the velocity 
of the frame we readily arrive at 

d 2 H dH 

(Hr + ifn - iu)H + (D r + iDi)—^ + v- {f3 r + if3i)\H\ 2 H = 

dy dy 



Let us now assume that H has a local form like He ^ . With this consideration the above 
equation changes to 

,2 ; 



[D r + Di 
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+ (fr + ifii - iu)He— + (f3 r + i^)H A e— = (1) 



Since we are interested in a form of H which is equal to ye ^ we can safely drop the 
nonlinear term in H. This approximation will be automatically justified latter. A point 

to be noted here is that, the linear solution is everywhere valid within the envelop e 26 
Since, near origin y is small H is small and if the envelop is not too wide, away from the 
origin H is also small and the nonlinear term do not count. 

The Equations coming out by equating real and imaginary parts will now look like 



D r 



d 2 H 2ydH y 2 H H 
'dy 2 ~ T~dy~ + ~W ~ J 



+ v 



dH Hy 
dy b 



+ jj, r H = 



and 



Di 



d 2 H 2ydH y 2 H H 



+ v 



dH Hy 

dy b 



(m-u)H = o 



They look almost the same except for the last terms. Now, considering one that comes 
from the real parts and rearranging the terms we get 

„ d 2 H 2D r . vb .dH D rt vb . 2tt b 2 v 2 tt , D r . TT n 

D ^W-^ (y -2iT r ^ + ^ y -2ij: )H -w; H+ ^-T )H = 



(2) 

The above form of the equation suggests a natural shift of the origin by an amount 
The shift in the origin will be more with faster moving frames and with a wider spread of 
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the localized solution. Let us effect this shift by going to z — y — and after effecting 

2 

that to look at around this new origin one can easily drop the term containing j~ with a 
promise that we will consider only linear terms in the solution of this equation. So, the 
region of interest has a spread of b on the far side of this shifted origin because of the fact 
that the nonlinear term in H is more relevant towards the inner side. Now, the equation 
takes the form 

n d 2 H 2D r zdH . D r b 2 v\ TT 
dz 2 b dz VPr b AD r J 



Now, the last transformation required is going to a larger scale defined by z = ^= 
which will produce the final form of the equation as 



d 2 H n dH b . D T b 2 v 2 STT n 
— — ^- - 2z— + —(ur - -rrr)H = 

dz 2 dz Dr b AD/ 



(3) 



The linear Eq.3 admits solutions which are Hermite polynomials and a solution of 
order zero is obtained when (fj, r — — is equal to zero. When fi r — ^- — = 2 -y L it 
admits a solution same as Hermite polynomial of order unity. We would not consider the 
solutions of higher order than unity since we have ignored the z 2 term to arrive at Eq.3 
but this will suffice to serve our purpose. Now, consider the steady state solution to Eq.3 
in the form of Hermite polynomial of order unity. This profile will be uniformly moving 
toward right starting from a point separated by right of the stationary origin at time 
t — 0. what about the same thing with v = —vl There will be a similar local structure 
moving toward left starting from the left of the stationary origin at a distance — at 
£ = 0. Thus, we can actually think about separating the space simultaneously into two 
moving frames on the two sides of the stationary core while on one a pulse moves toward 
right and on the other it moves toward left. 

Let us think about the effect of the shift in origin by the amount from the center of 
=£. 

the envelop given by e ^ in the case of the right-moving solution. In this case, within the 
envelop the right hand part of the Hermite polynomial of order unity will be suppressed 
and the left hand part will be most prominent. So, There will be a trough moving towards 
right hand side when by the same logic, due to the shift in origin within the left-moving 
envelop, a peak will be moving towards the laft hand side. It is interesting to see that 
for a narrow envelop (b small) the suppression of one half of the spatial solution is more 
prominent and the antisymmetric nature of the solution will be more revealed and at 
the same time a narrow envelop makes the nonlinear effect even more negligible. Now, if 
we remember the fact that these structures are oscillating with frequency u, we readily 
realize that this scenario will be changing with time. So, a peak will follow a trough on 
the right hand side while a trough will follow a peak on the left after a time 

Considering the equation obtained from imaginary part along with that from real part 
requires simultaneous validity of 

A b 2 v 2 , 2 A 
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So, the solution now can be obtained when 

b b 3 v 2 b b 3 v 2 

meaning that, on the straight line in fi r — \ii plane for given other parameters, u> being the 
frequency of oscillation of moving local structures can adjust itself to lock to the solution 
of order unity. 

In the situation when D r and are negative, one has to add positive quartic terms 
to the amplitude equation to make the system saturate. Thus, in such a situation the 
system is sub critical. Since now changes sign the right moving pulse will have its 
origin on the left of the zero of y and vice versa. In this situation the assumption of two 
independent frames on either sides breaks down and as a result such phenomenon are 
strictly supercritical. We would like to mention that super criticality has also been taken 
as a condition in numerical simulations to show generation of ID spirals on the basis 
of multi stability pQ. As the present theory shows, one can really expect such moving 
local structures can also be seen in a pure Turing environment near instability boundary 
where the same form of amplitude equation applies but with real coefficients. The only 
difference will be that those moving localized structures will be non oscillatory and thus 
concentration peaks will travel toward right while troughs toward left. Experiments, in 
this direction, to capture such non oscillatory peaks and troughs moving on the two sides 
of a core can show the validity of this theory and can verify if the non variational effects 
are essential or not for such pulse generation. 

In the conclusion we would like to say that we have put forward an alternative theory 
for antisymmetric pulse propagation or ID spirals which has nothing to do with multi 
stability of states. The consideration of non variational effect comes out to be nonessential 
in view of the fact that such pulses can arise in a pure Turing type amplitude equation 
which has a potential form. There are certainly inherent mechanism other than multi 
stability of phases which is also capable of showing the phenomenon of ID spiral. 
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